Here, I’m trying to get a better sense of how well the SODA temp data represents the SEUS. We found a significant decrease in temperature during the study period, which is not what people who are familiar with conditions in this area say has happened. Most likely, the SODA model doesn’t perform well in shallow near shore areas where survey is performed in southeast US. Therefore, we compared performance to OISST temp. Overall pattern is similar, and not representative of what’s actually happening in this survey area. Therefore, we decided to go with in situ buoy temperature values.

library(ncdf4)
library(raster)
Loading required package: sp

Attaching package: ‘raster’

The following object is masked from ‘package:dplyr’:

    select

The following object is masked from ‘package:tidyr’:

    extract

The following object is masked from ‘package:lme4’:

    getData

The following object is masked from ‘package:data.table’:

    shift
library(data.table)

Southeast US (1990–2014, n = 25)

Years downloaded one by one on May 15, 2020 https://psl.noaa.gov/cgi-bin/DataAccess.pl?DB_dataset=NOAA+High-resolution+Blended+Analysis&DB_variable=Sea+Surface+Temperature&DB_statistic=Mean&DB_tid=86107&DB_did=132&DB_vid=2423

#all file names
filenames <- list.files(here::here("Data","Temperature","OISST_NOAA_TEMP_NETCDF"), full.names = T)

OISST.sst <- lapply(filenames, brick)

OISST.sst
[[1]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1990.nc 
names      : X1990.01.01, X1990.01.02, X1990.01.03, X1990.01.04, X1990.01.05, X1990.01.06, X1990.01.07, X1990.01.08, X1990.01.09, X1990.01.10, X1990.01.11, X1990.01.12, X1990.01.13, X1990.01.14, X1990.01.15, ... 
Date       : 1990-01-01, 1990-12-31 (min, max)
varname    : sst 


[[2]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1991.nc 
names      : X1991.01.01, X1991.01.02, X1991.01.03, X1991.01.04, X1991.01.05, X1991.01.06, X1991.01.07, X1991.01.08, X1991.01.09, X1991.01.10, X1991.01.11, X1991.01.12, X1991.01.13, X1991.01.14, X1991.01.15, ... 
Date       : 1991-01-01, 1991-12-31 (min, max)
varname    : sst 


[[3]]
class      : RasterBrick 
dimensions : 32, 28, 896, 366  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1992.nc 
names      : X1992.01.01, X1992.01.02, X1992.01.03, X1992.01.04, X1992.01.05, X1992.01.06, X1992.01.07, X1992.01.08, X1992.01.09, X1992.01.10, X1992.01.11, X1992.01.12, X1992.01.13, X1992.01.14, X1992.01.15, ... 
Date       : 1992-01-01, 1992-12-31 (min, max)
varname    : sst 


[[4]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1993.nc 
names      : X1993.01.01, X1993.01.02, X1993.01.03, X1993.01.04, X1993.01.05, X1993.01.06, X1993.01.07, X1993.01.08, X1993.01.09, X1993.01.10, X1993.01.11, X1993.01.12, X1993.01.13, X1993.01.14, X1993.01.15, ... 
Date       : 1993-01-01, 1993-12-31 (min, max)
varname    : sst 


[[5]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1994.nc 
names      : X1994.01.01, X1994.01.02, X1994.01.03, X1994.01.04, X1994.01.05, X1994.01.06, X1994.01.07, X1994.01.08, X1994.01.09, X1994.01.10, X1994.01.11, X1994.01.12, X1994.01.13, X1994.01.14, X1994.01.15, ... 
Date       : 1994-01-01, 1994-12-31 (min, max)
varname    : sst 


[[6]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1995.nc 
names      : X1995.01.01, X1995.01.02, X1995.01.03, X1995.01.04, X1995.01.05, X1995.01.06, X1995.01.07, X1995.01.08, X1995.01.09, X1995.01.10, X1995.01.11, X1995.01.12, X1995.01.13, X1995.01.14, X1995.01.15, ... 
Date       : 1995-01-01, 1995-12-31 (min, max)
varname    : sst 


[[7]]
class      : RasterBrick 
dimensions : 32, 28, 896, 366  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1996.nc 
names      : X1996.01.01, X1996.01.02, X1996.01.03, X1996.01.04, X1996.01.05, X1996.01.06, X1996.01.07, X1996.01.08, X1996.01.09, X1996.01.10, X1996.01.11, X1996.01.12, X1996.01.13, X1996.01.14, X1996.01.15, ... 
Date       : 1996-01-01, 1996-12-31 (min, max)
varname    : sst 


[[8]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1997.nc 
names      : X1997.01.01, X1997.01.02, X1997.01.03, X1997.01.04, X1997.01.05, X1997.01.06, X1997.01.07, X1997.01.08, X1997.01.09, X1997.01.10, X1997.01.11, X1997.01.12, X1997.01.13, X1997.01.14, X1997.01.15, ... 
Date       : 1997-01-01, 1997-12-31 (min, max)
varname    : sst 


[[9]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1998.nc 
names      : X1998.01.01, X1998.01.02, X1998.01.03, X1998.01.04, X1998.01.05, X1998.01.06, X1998.01.07, X1998.01.08, X1998.01.09, X1998.01.10, X1998.01.11, X1998.01.12, X1998.01.13, X1998.01.14, X1998.01.15, ... 
Date       : 1998-01-01, 1998-12-31 (min, max)
varname    : sst 


[[10]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1999.nc 
names      : X1999.01.01, X1999.01.02, X1999.01.03, X1999.01.04, X1999.01.05, X1999.01.06, X1999.01.07, X1999.01.08, X1999.01.09, X1999.01.10, X1999.01.11, X1999.01.12, X1999.01.13, X1999.01.14, X1999.01.15, ... 
Date       : 1999-01-01, 1999-12-31 (min, max)
varname    : sst 


[[11]]
class      : RasterBrick 
dimensions : 32, 28, 896, 366  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2000.nc 
names      : X2000.01.01, X2000.01.02, X2000.01.03, X2000.01.04, X2000.01.05, X2000.01.06, X2000.01.07, X2000.01.08, X2000.01.09, X2000.01.10, X2000.01.11, X2000.01.12, X2000.01.13, X2000.01.14, X2000.01.15, ... 
Date       : 2000-01-01, 2000-12-31 (min, max)
varname    : sst 


[[12]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2001.nc 
names      : X2001.01.01, X2001.01.02, X2001.01.03, X2001.01.04, X2001.01.05, X2001.01.06, X2001.01.07, X2001.01.08, X2001.01.09, X2001.01.10, X2001.01.11, X2001.01.12, X2001.01.13, X2001.01.14, X2001.01.15, ... 
Date       : 2001-01-01, 2001-12-31 (min, max)
varname    : sst 


[[13]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2002.nc 
names      : X2002.01.01, X2002.01.02, X2002.01.03, X2002.01.04, X2002.01.05, X2002.01.06, X2002.01.07, X2002.01.08, X2002.01.09, X2002.01.10, X2002.01.11, X2002.01.12, X2002.01.13, X2002.01.14, X2002.01.15, ... 
Date       : 2002-01-01, 2002-12-31 (min, max)
varname    : sst 


[[14]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2003.nc 
names      : X2003.01.01, X2003.01.02, X2003.01.03, X2003.01.04, X2003.01.05, X2003.01.06, X2003.01.07, X2003.01.08, X2003.01.09, X2003.01.10, X2003.01.11, X2003.01.12, X2003.01.13, X2003.01.14, X2003.01.15, ... 
Date       : 2003-01-01, 2003-12-31 (min, max)
varname    : sst 


[[15]]
class      : RasterBrick 
dimensions : 32, 28, 896, 366  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2004.nc 
names      : X2004.01.01, X2004.01.02, X2004.01.03, X2004.01.04, X2004.01.05, X2004.01.06, X2004.01.07, X2004.01.08, X2004.01.09, X2004.01.10, X2004.01.11, X2004.01.12, X2004.01.13, X2004.01.14, X2004.01.15, ... 
Date       : 2004-01-01, 2004-12-31 (min, max)
varname    : sst 


[[16]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2005.nc 
names      : X2005.01.01, X2005.01.02, X2005.01.03, X2005.01.04, X2005.01.05, X2005.01.06, X2005.01.07, X2005.01.08, X2005.01.09, X2005.01.10, X2005.01.11, X2005.01.12, X2005.01.13, X2005.01.14, X2005.01.15, ... 
Date       : 2005-01-01, 2005-12-31 (min, max)
varname    : sst 


[[17]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2006.nc 
names      : X2006.01.01, X2006.01.02, X2006.01.03, X2006.01.04, X2006.01.05, X2006.01.06, X2006.01.07, X2006.01.08, X2006.01.09, X2006.01.10, X2006.01.11, X2006.01.12, X2006.01.13, X2006.01.14, X2006.01.15, ... 
Date       : 2006-01-01, 2006-12-31 (min, max)
varname    : sst 


[[18]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2007.nc 
names      : X2007.01.01, X2007.01.02, X2007.01.03, X2007.01.04, X2007.01.05, X2007.01.06, X2007.01.07, X2007.01.08, X2007.01.09, X2007.01.10, X2007.01.11, X2007.01.12, X2007.01.13, X2007.01.14, X2007.01.15, ... 
Date       : 2007-01-01, 2007-12-31 (min, max)
varname    : sst 


[[19]]
class      : RasterBrick 
dimensions : 32, 28, 896, 366  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2008.nc 
names      : X2008.01.01, X2008.01.02, X2008.01.03, X2008.01.04, X2008.01.05, X2008.01.06, X2008.01.07, X2008.01.08, X2008.01.09, X2008.01.10, X2008.01.11, X2008.01.12, X2008.01.13, X2008.01.14, X2008.01.15, ... 
Date       : 2008-01-01, 2008-12-31 (min, max)
varname    : sst 


[[20]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2009.nc 
names      : X2009.01.01, X2009.01.02, X2009.01.03, X2009.01.04, X2009.01.05, X2009.01.06, X2009.01.07, X2009.01.08, X2009.01.09, X2009.01.10, X2009.01.11, X2009.01.12, X2009.01.13, X2009.01.14, X2009.01.15, ... 
Date       : 2009-01-01, 2009-12-31 (min, max)
varname    : sst 


[[21]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2010.nc 
names      : X2010.01.01, X2010.01.02, X2010.01.03, X2010.01.04, X2010.01.05, X2010.01.06, X2010.01.07, X2010.01.08, X2010.01.09, X2010.01.10, X2010.01.11, X2010.01.12, X2010.01.13, X2010.01.14, X2010.01.15, ... 
Date       : 2010-01-01, 2010-12-31 (min, max)
varname    : sst 


[[22]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2011.nc 
names      : X2011.01.01, X2011.01.02, X2011.01.03, X2011.01.04, X2011.01.05, X2011.01.06, X2011.01.07, X2011.01.08, X2011.01.09, X2011.01.10, X2011.01.11, X2011.01.12, X2011.01.13, X2011.01.14, X2011.01.15, ... 
Date       : 2011-01-01, 2011-12-31 (min, max)
varname    : sst 


[[23]]
class      : RasterBrick 
dimensions : 32, 28, 896, 366  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2012.nc 
names      : X2012.01.01, X2012.01.02, X2012.01.03, X2012.01.04, X2012.01.05, X2012.01.06, X2012.01.07, X2012.01.08, X2012.01.09, X2012.01.10, X2012.01.11, X2012.01.12, X2012.01.13, X2012.01.14, X2012.01.15, ... 
Date       : 2012-01-01, 2012-12-31 (min, max)
varname    : sst 


[[24]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2013.nc 
names      : X2013.01.01, X2013.01.02, X2013.01.03, X2013.01.04, X2013.01.05, X2013.01.06, X2013.01.07, X2013.01.08, X2013.01.09, X2013.01.10, X2013.01.11, X2013.01.12, X2013.01.13, X2013.01.14, X2013.01.15, ... 
Date       : 2013-01-01, 2013-12-31 (min, max)
varname    : sst 


[[25]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2014.nc 
names      : X2014.01.01, X2014.01.02, X2014.01.03, X2014.01.04, X2014.01.05, X2014.01.06, X2014.01.07, X2014.01.08, X2014.01.09, X2014.01.10, X2014.01.11, X2014.01.12, X2014.01.13, X2014.01.14, X2014.01.15, ... 
Date       : 2014-01-01, 2014-12-31 (min, max)
varname    : sst 


[[26]]
class      : RasterBrick 
dimensions : 32, 28, 896, 365  (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25  (x, y)
extent     : 278, 285, 28, 36  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
source     : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2015.nc 
names      : X2015.01.01, X2015.01.02, X2015.01.03, X2015.01.04, X2015.01.05, X2015.01.06, X2015.01.07, X2015.01.08, X2015.01.09, X2015.01.10, X2015.01.11, X2015.01.12, X2015.01.13, X2015.01.14, X2015.01.15, ... 
Date       : 2015-01-01, 2015-12-31 (min, max)
varname    : sst 
#currently, this is in 0-360, needs to be -180-180, so I have to subtract 360Ëš
OISST.sst_adj <- lapply(OISST.sst, rotate)
this does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this function

Extract specific GPS values for SE US sampled during during study period

load(here::here("Data","Spp_observations","trawl_coordinates_1996-2000.RData"))

seus_points <- gpspointsonly_quarterdegreeround[region == "SCDNR_SEUS"]

coordinates(seus_points) = ~ lon + lat


temp_values_1990 <- as.data.table(extract(OISST.sst_adj[[1]], seus_points))
temp_values_1991 <- as.data.table(extract(OISST.sst_adj[[2]], seus_points))
temp_values_1992 <- as.data.table(extract(OISST.sst_adj[[3]], seus_points))
temp_values_1993 <- as.data.table(extract(OISST.sst_adj[[4]], seus_points))
temp_values_1994 <- as.data.table(extract(OISST.sst_adj[[5]], seus_points))
temp_values_1995 <- as.data.table(extract(OISST.sst_adj[[6]], seus_points))
temp_values_1996 <- as.data.table(extract(OISST.sst_adj[[7]], seus_points))
temp_values_1997 <- as.data.table(extract(OISST.sst_adj[[8]], seus_points))
temp_values_1998 <- as.data.table(extract(OISST.sst_adj[[9]], seus_points))
temp_values_1999 <- as.data.table(extract(OISST.sst_adj[[10]], seus_points))
temp_values_2000 <- as.data.table(extract(OISST.sst_adj[[11]], seus_points))
temp_values_2001 <- as.data.table(extract(OISST.sst_adj[[12]], seus_points))
temp_values_2002 <- as.data.table(extract(OISST.sst_adj[[13]], seus_points))
temp_values_2003 <- as.data.table(extract(OISST.sst_adj[[14]], seus_points))
temp_values_2004 <- as.data.table(extract(OISST.sst_adj[[15]], seus_points))
temp_values_2005 <- as.data.table(extract(OISST.sst_adj[[16]], seus_points))
temp_values_2006 <- as.data.table(extract(OISST.sst_adj[[17]], seus_points))
temp_values_2007 <- as.data.table(extract(OISST.sst_adj[[18]], seus_points))
temp_values_2008 <- as.data.table(extract(OISST.sst_adj[[19]], seus_points))
temp_values_2009 <- as.data.table(extract(OISST.sst_adj[[20]], seus_points))
temp_values_2010 <- as.data.table(extract(OISST.sst_adj[[21]], seus_points))
temp_values_2011 <- as.data.table(extract(OISST.sst_adj[[22]], seus_points))
temp_values_2012 <- as.data.table(extract(OISST.sst_adj[[23]], seus_points))
temp_values_2013 <- as.data.table(extract(OISST.sst_adj[[24]], seus_points))
temp_values_2014 <- as.data.table(extract(OISST.sst_adj[[25]], seus_points))
temp_values_2015 <- as.data.table(extract(OISST.sst_adj[[26]], seus_points))

plot(OISST.sst_adj[[1]][[1]])
points(seus_points)


object_names <- list(temp_values_1990, temp_values_1991, temp_values_1992, temp_values_1993, temp_values_1994, temp_values_1995, temp_values_1996, temp_values_1997, temp_values_1998, temp_values_1999, temp_values_2000, temp_values_2001, temp_values_2002, temp_values_2003, temp_values_2004, temp_values_2005, temp_values_2006, temp_values_2007, temp_values_2008, temp_values_2009, temp_values_2010, temp_values_2011, temp_values_2012, temp_values_2013, temp_values_2014, temp_values_2015)

allvalues_wide <- do.call("cbind", object_names)

allvalues_wide_latlon <- cbind(seus_points@coords, allvalues_wide)
colnames(allvalues_wide_latlon)
   [1] "lon"         "lat"         "X1990.01.01" "X1990.01.02" "X1990.01.03" "X1990.01.04" "X1990.01.05"
   [8] "X1990.01.06" "X1990.01.07" "X1990.01.08" "X1990.01.09" "X1990.01.10" "X1990.01.11" "X1990.01.12"
  [15] "X1990.01.13" "X1990.01.14" "X1990.01.15" "X1990.01.16" "X1990.01.17" "X1990.01.18" "X1990.01.19"
  [22] "X1990.01.20" "X1990.01.21" "X1990.01.22" "X1990.01.23" "X1990.01.24" "X1990.01.25" "X1990.01.26"
  [29] "X1990.01.27" "X1990.01.28" "X1990.01.29" "X1990.01.30" "X1990.01.31" "X1990.02.01" "X1990.02.02"
  [36] "X1990.02.03" "X1990.02.04" "X1990.02.05" "X1990.02.06" "X1990.02.07" "X1990.02.08" "X1990.02.09"
  [43] "X1990.02.10" "X1990.02.11" "X1990.02.12" "X1990.02.13" "X1990.02.14" "X1990.02.15" "X1990.02.16"
  [50] "X1990.02.17" "X1990.02.18" "X1990.02.19" "X1990.02.20" "X1990.02.21" "X1990.02.22" "X1990.02.23"
  [57] "X1990.02.24" "X1990.02.25" "X1990.02.26" "X1990.02.27" "X1990.02.28" "X1990.03.01" "X1990.03.02"
  [64] "X1990.03.03" "X1990.03.04" "X1990.03.05" "X1990.03.06" "X1990.03.07" "X1990.03.08" "X1990.03.09"
  [71] "X1990.03.10" "X1990.03.11" "X1990.03.12" "X1990.03.13" "X1990.03.14" "X1990.03.15" "X1990.03.16"
  [78] "X1990.03.17" "X1990.03.18" "X1990.03.19" "X1990.03.20" "X1990.03.21" "X1990.03.22" "X1990.03.23"
  [85] "X1990.03.24" "X1990.03.25" "X1990.03.26" "X1990.03.27" "X1990.03.28" "X1990.03.29" "X1990.03.30"
  [92] "X1990.03.31" "X1990.04.01" "X1990.04.02" "X1990.04.03" "X1990.04.04" "X1990.04.05" "X1990.04.06"
  [99] "X1990.04.07" "X1990.04.08" "X1990.04.09" "X1990.04.10" "X1990.04.11" "X1990.04.12" "X1990.04.13"
 [106] "X1990.04.14" "X1990.04.15" "X1990.04.16" "X1990.04.17" "X1990.04.18" "X1990.04.19" "X1990.04.20"
 [113] "X1990.04.21" "X1990.04.22" "X1990.04.23" "X1990.04.24" "X1990.04.25" "X1990.04.26" "X1990.04.27"
 [120] "X1990.04.28" "X1990.04.29" "X1990.04.30" "X1990.05.01" "X1990.05.02" "X1990.05.03" "X1990.05.04"
 [127] "X1990.05.05" "X1990.05.06" "X1990.05.07" "X1990.05.08" "X1990.05.09" "X1990.05.10" "X1990.05.11"
 [134] "X1990.05.12" "X1990.05.13" "X1990.05.14" "X1990.05.15" "X1990.05.16" "X1990.05.17" "X1990.05.18"
 [141] "X1990.05.19" "X1990.05.20" "X1990.05.21" "X1990.05.22" "X1990.05.23" "X1990.05.24" "X1990.05.25"
 [148] "X1990.05.26" "X1990.05.27" "X1990.05.28" "X1990.05.29" "X1990.05.30" "X1990.05.31" "X1990.06.01"
 [155] "X1990.06.02" "X1990.06.03" "X1990.06.04" "X1990.06.05" "X1990.06.06" "X1990.06.07" "X1990.06.08"
 [162] "X1990.06.09" "X1990.06.10" "X1990.06.11" "X1990.06.12" "X1990.06.13" "X1990.06.14" "X1990.06.15"
 [169] "X1990.06.16" "X1990.06.17" "X1990.06.18" "X1990.06.19" "X1990.06.20" "X1990.06.21" "X1990.06.22"
 [176] "X1990.06.23" "X1990.06.24" "X1990.06.25" "X1990.06.26" "X1990.06.27" "X1990.06.28" "X1990.06.29"
 [183] "X1990.06.30" "X1990.07.01" "X1990.07.02" "X1990.07.03" "X1990.07.04" "X1990.07.05" "X1990.07.06"
 [190] "X1990.07.07" "X1990.07.08" "X1990.07.09" "X1990.07.10" "X1990.07.11" "X1990.07.12" "X1990.07.13"
 [197] "X1990.07.14" "X1990.07.15" "X1990.07.16" "X1990.07.17" "X1990.07.18" "X1990.07.19" "X1990.07.20"
 [204] "X1990.07.21" "X1990.07.22" "X1990.07.23" "X1990.07.24" "X1990.07.25" "X1990.07.26" "X1990.07.27"
 [211] "X1990.07.28" "X1990.07.29" "X1990.07.30" "X1990.07.31" "X1990.08.01" "X1990.08.02" "X1990.08.03"
 [218] "X1990.08.04" "X1990.08.05" "X1990.08.06" "X1990.08.07" "X1990.08.08" "X1990.08.09" "X1990.08.10"
 [225] "X1990.08.11" "X1990.08.12" "X1990.08.13" "X1990.08.14" "X1990.08.15" "X1990.08.16" "X1990.08.17"
 [232] "X1990.08.18" "X1990.08.19" "X1990.08.20" "X1990.08.21" "X1990.08.22" "X1990.08.23" "X1990.08.24"
 [239] "X1990.08.25" "X1990.08.26" "X1990.08.27" "X1990.08.28" "X1990.08.29" "X1990.08.30" "X1990.08.31"
 [246] "X1990.09.01" "X1990.09.02" "X1990.09.03" "X1990.09.04" "X1990.09.05" "X1990.09.06" "X1990.09.07"
 [253] "X1990.09.08" "X1990.09.09" "X1990.09.10" "X1990.09.11" "X1990.09.12" "X1990.09.13" "X1990.09.14"
 [260] "X1990.09.15" "X1990.09.16" "X1990.09.17" "X1990.09.18" "X1990.09.19" "X1990.09.20" "X1990.09.21"
 [267] "X1990.09.22" "X1990.09.23" "X1990.09.24" "X1990.09.25" "X1990.09.26" "X1990.09.27" "X1990.09.28"
 [274] "X1990.09.29" "X1990.09.30" "X1990.10.01" "X1990.10.02" "X1990.10.03" "X1990.10.04" "X1990.10.05"
 [281] "X1990.10.06" "X1990.10.07" "X1990.10.08" "X1990.10.09" "X1990.10.10" "X1990.10.11" "X1990.10.12"
 [288] "X1990.10.13" "X1990.10.14" "X1990.10.15" "X1990.10.16" "X1990.10.17" "X1990.10.18" "X1990.10.19"
 [295] "X1990.10.20" "X1990.10.21" "X1990.10.22" "X1990.10.23" "X1990.10.24" "X1990.10.25" "X1990.10.26"
 [302] "X1990.10.27" "X1990.10.28" "X1990.10.29" "X1990.10.30" "X1990.10.31" "X1990.11.01" "X1990.11.02"
 [309] "X1990.11.03" "X1990.11.04" "X1990.11.05" "X1990.11.06" "X1990.11.07" "X1990.11.08" "X1990.11.09"
 [316] "X1990.11.10" "X1990.11.11" "X1990.11.12" "X1990.11.13" "X1990.11.14" "X1990.11.15" "X1990.11.16"
 [323] "X1990.11.17" "X1990.11.18" "X1990.11.19" "X1990.11.20" "X1990.11.21" "X1990.11.22" "X1990.11.23"
 [330] "X1990.11.24" "X1990.11.25" "X1990.11.26" "X1990.11.27" "X1990.11.28" "X1990.11.29" "X1990.11.30"
 [337] "X1990.12.01" "X1990.12.02" "X1990.12.03" "X1990.12.04" "X1990.12.05" "X1990.12.06" "X1990.12.07"
 [344] "X1990.12.08" "X1990.12.09" "X1990.12.10" "X1990.12.11" "X1990.12.12" "X1990.12.13" "X1990.12.14"
 [351] "X1990.12.15" "X1990.12.16" "X1990.12.17" "X1990.12.18" "X1990.12.19" "X1990.12.20" "X1990.12.21"
 [358] "X1990.12.22" "X1990.12.23" "X1990.12.24" "X1990.12.25" "X1990.12.26" "X1990.12.27" "X1990.12.28"
 [365] "X1990.12.29" "X1990.12.30" "X1990.12.31" "X1991.01.01" "X1991.01.02" "X1991.01.03" "X1991.01.04"
 [372] "X1991.01.05" "X1991.01.06" "X1991.01.07" "X1991.01.08" "X1991.01.09" "X1991.01.10" "X1991.01.11"
 [379] "X1991.01.12" "X1991.01.13" "X1991.01.14" "X1991.01.15" "X1991.01.16" "X1991.01.17" "X1991.01.18"
 [386] "X1991.01.19" "X1991.01.20" "X1991.01.21" "X1991.01.22" "X1991.01.23" "X1991.01.24" "X1991.01.25"
 [393] "X1991.01.26" "X1991.01.27" "X1991.01.28" "X1991.01.29" "X1991.01.30" "X1991.01.31" "X1991.02.01"
 [400] "X1991.02.02" "X1991.02.03" "X1991.02.04" "X1991.02.05" "X1991.02.06" "X1991.02.07" "X1991.02.08"
 [407] "X1991.02.09" "X1991.02.10" "X1991.02.11" "X1991.02.12" "X1991.02.13" "X1991.02.14" "X1991.02.15"
 [414] "X1991.02.16" "X1991.02.17" "X1991.02.18" "X1991.02.19" "X1991.02.20" "X1991.02.21" "X1991.02.22"
 [421] "X1991.02.23" "X1991.02.24" "X1991.02.25" "X1991.02.26" "X1991.02.27" "X1991.02.28" "X1991.03.01"
 [428] "X1991.03.02" "X1991.03.03" "X1991.03.04" "X1991.03.05" "X1991.03.06" "X1991.03.07" "X1991.03.08"
 [435] "X1991.03.09" "X1991.03.10" "X1991.03.11" "X1991.03.12" "X1991.03.13" "X1991.03.14" "X1991.03.15"
 [442] "X1991.03.16" "X1991.03.17" "X1991.03.18" "X1991.03.19" "X1991.03.20" "X1991.03.21" "X1991.03.22"
 [449] "X1991.03.23" "X1991.03.24" "X1991.03.25" "X1991.03.26" "X1991.03.27" "X1991.03.28" "X1991.03.29"
 [456] "X1991.03.30" "X1991.03.31" "X1991.04.01" "X1991.04.02" "X1991.04.03" "X1991.04.04" "X1991.04.05"
 [463] "X1991.04.06" "X1991.04.07" "X1991.04.08" "X1991.04.09" "X1991.04.10" "X1991.04.11" "X1991.04.12"
 [470] "X1991.04.13" "X1991.04.14" "X1991.04.15" "X1991.04.16" "X1991.04.17" "X1991.04.18" "X1991.04.19"
 [477] "X1991.04.20" "X1991.04.21" "X1991.04.22" "X1991.04.23" "X1991.04.24" "X1991.04.25" "X1991.04.26"
 [484] "X1991.04.27" "X1991.04.28" "X1991.04.29" "X1991.04.30" "X1991.05.01" "X1991.05.02" "X1991.05.03"
 [491] "X1991.05.04" "X1991.05.05" "X1991.05.06" "X1991.05.07" "X1991.05.08" "X1991.05.09" "X1991.05.10"
 [498] "X1991.05.11" "X1991.05.12" "X1991.05.13" "X1991.05.14" "X1991.05.15" "X1991.05.16" "X1991.05.17"
 [505] "X1991.05.18" "X1991.05.19" "X1991.05.20" "X1991.05.21" "X1991.05.22" "X1991.05.23" "X1991.05.24"
 [512] "X1991.05.25" "X1991.05.26" "X1991.05.27" "X1991.05.28" "X1991.05.29" "X1991.05.30" "X1991.05.31"
 [519] "X1991.06.01" "X1991.06.02" "X1991.06.03" "X1991.06.04" "X1991.06.05" "X1991.06.06" "X1991.06.07"
 [526] "X1991.06.08" "X1991.06.09" "X1991.06.10" "X1991.06.11" "X1991.06.12" "X1991.06.13" "X1991.06.14"
 [533] "X1991.06.15" "X1991.06.16" "X1991.06.17" "X1991.06.18" "X1991.06.19" "X1991.06.20" "X1991.06.21"
 [540] "X1991.06.22" "X1991.06.23" "X1991.06.24" "X1991.06.25" "X1991.06.26" "X1991.06.27" "X1991.06.28"
 [547] "X1991.06.29" "X1991.06.30" "X1991.07.01" "X1991.07.02" "X1991.07.03" "X1991.07.04" "X1991.07.05"
 [554] "X1991.07.06" "X1991.07.07" "X1991.07.08" "X1991.07.09" "X1991.07.10" "X1991.07.11" "X1991.07.12"
 [561] "X1991.07.13" "X1991.07.14" "X1991.07.15" "X1991.07.16" "X1991.07.17" "X1991.07.18" "X1991.07.19"
 [568] "X1991.07.20" "X1991.07.21" "X1991.07.22" "X1991.07.23" "X1991.07.24" "X1991.07.25" "X1991.07.26"
 [575] "X1991.07.27" "X1991.07.28" "X1991.07.29" "X1991.07.30" "X1991.07.31" "X1991.08.01" "X1991.08.02"
 [582] "X1991.08.03" "X1991.08.04" "X1991.08.05" "X1991.08.06" "X1991.08.07" "X1991.08.08" "X1991.08.09"
 [589] "X1991.08.10" "X1991.08.11" "X1991.08.12" "X1991.08.13" "X1991.08.14" "X1991.08.15" "X1991.08.16"
 [596] "X1991.08.17" "X1991.08.18" "X1991.08.19" "X1991.08.20" "X1991.08.21" "X1991.08.22" "X1991.08.23"
 [603] "X1991.08.24" "X1991.08.25" "X1991.08.26" "X1991.08.27" "X1991.08.28" "X1991.08.29" "X1991.08.30"
 [610] "X1991.08.31" "X1991.09.01" "X1991.09.02" "X1991.09.03" "X1991.09.04" "X1991.09.05" "X1991.09.06"
 [617] "X1991.09.07" "X1991.09.08" "X1991.09.09" "X1991.09.10" "X1991.09.11" "X1991.09.12" "X1991.09.13"
 [624] "X1991.09.14" "X1991.09.15" "X1991.09.16" "X1991.09.17" "X1991.09.18" "X1991.09.19" "X1991.09.20"
 [631] "X1991.09.21" "X1991.09.22" "X1991.09.23" "X1991.09.24" "X1991.09.25" "X1991.09.26" "X1991.09.27"
 [638] "X1991.09.28" "X1991.09.29" "X1991.09.30" "X1991.10.01" "X1991.10.02" "X1991.10.03" "X1991.10.04"
 [645] "X1991.10.05" "X1991.10.06" "X1991.10.07" "X1991.10.08" "X1991.10.09" "X1991.10.10" "X1991.10.11"
 [652] "X1991.10.12" "X1991.10.13" "X1991.10.14" "X1991.10.15" "X1991.10.16" "X1991.10.17" "X1991.10.18"
 [659] "X1991.10.19" "X1991.10.20" "X1991.10.21" "X1991.10.22" "X1991.10.23" "X1991.10.24" "X1991.10.25"
 [666] "X1991.10.26" "X1991.10.27" "X1991.10.28" "X1991.10.29" "X1991.10.30" "X1991.10.31" "X1991.11.01"
 [673] "X1991.11.02" "X1991.11.03" "X1991.11.04" "X1991.11.05" "X1991.11.06" "X1991.11.07" "X1991.11.08"
 [680] "X1991.11.09" "X1991.11.10" "X1991.11.11" "X1991.11.12" "X1991.11.13" "X1991.11.14" "X1991.11.15"
 [687] "X1991.11.16" "X1991.11.17" "X1991.11.18" "X1991.11.19" "X1991.11.20" "X1991.11.21" "X1991.11.22"
 [694] "X1991.11.23" "X1991.11.24" "X1991.11.25" "X1991.11.26" "X1991.11.27" "X1991.11.28" "X1991.11.29"
 [701] "X1991.11.30" "X1991.12.01" "X1991.12.02" "X1991.12.03" "X1991.12.04" "X1991.12.05" "X1991.12.06"
 [708] "X1991.12.07" "X1991.12.08" "X1991.12.09" "X1991.12.10" "X1991.12.11" "X1991.12.12" "X1991.12.13"
 [715] "X1991.12.14" "X1991.12.15" "X1991.12.16" "X1991.12.17" "X1991.12.18" "X1991.12.19" "X1991.12.20"
 [722] "X1991.12.21" "X1991.12.22" "X1991.12.23" "X1991.12.24" "X1991.12.25" "X1991.12.26" "X1991.12.27"
 [729] "X1991.12.28" "X1991.12.29" "X1991.12.30" "X1991.12.31" "X1992.01.01" "X1992.01.02" "X1992.01.03"
 [736] "X1992.01.04" "X1992.01.05" "X1992.01.06" "X1992.01.07" "X1992.01.08" "X1992.01.09" "X1992.01.10"
 [743] "X1992.01.11" "X1992.01.12" "X1992.01.13" "X1992.01.14" "X1992.01.15" "X1992.01.16" "X1992.01.17"
 [750] "X1992.01.18" "X1992.01.19" "X1992.01.20" "X1992.01.21" "X1992.01.22" "X1992.01.23" "X1992.01.24"
 [757] "X1992.01.25" "X1992.01.26" "X1992.01.27" "X1992.01.28" "X1992.01.29" "X1992.01.30" "X1992.01.31"
 [764] "X1992.02.01" "X1992.02.02" "X1992.02.03" "X1992.02.04" "X1992.02.05" "X1992.02.06" "X1992.02.07"
 [771] "X1992.02.08" "X1992.02.09" "X1992.02.10" "X1992.02.11" "X1992.02.12" "X1992.02.13" "X1992.02.14"
 [778] "X1992.02.15" "X1992.02.16" "X1992.02.17" "X1992.02.18" "X1992.02.19" "X1992.02.20" "X1992.02.21"
 [785] "X1992.02.22" "X1992.02.23" "X1992.02.24" "X1992.02.25" "X1992.02.26" "X1992.02.27" "X1992.02.28"
 [792] "X1992.02.29" "X1992.03.01" "X1992.03.02" "X1992.03.03" "X1992.03.04" "X1992.03.05" "X1992.03.06"
 [799] "X1992.03.07" "X1992.03.08" "X1992.03.09" "X1992.03.10" "X1992.03.11" "X1992.03.12" "X1992.03.13"
 [806] "X1992.03.14" "X1992.03.15" "X1992.03.16" "X1992.03.17" "X1992.03.18" "X1992.03.19" "X1992.03.20"
 [813] "X1992.03.21" "X1992.03.22" "X1992.03.23" "X1992.03.24" "X1992.03.25" "X1992.03.26" "X1992.03.27"
 [820] "X1992.03.28" "X1992.03.29" "X1992.03.30" "X1992.03.31" "X1992.04.01" "X1992.04.02" "X1992.04.03"
 [827] "X1992.04.04" "X1992.04.05" "X1992.04.06" "X1992.04.07" "X1992.04.08" "X1992.04.09" "X1992.04.10"
 [834] "X1992.04.11" "X1992.04.12" "X1992.04.13" "X1992.04.14" "X1992.04.15" "X1992.04.16" "X1992.04.17"
 [841] "X1992.04.18" "X1992.04.19" "X1992.04.20" "X1992.04.21" "X1992.04.22" "X1992.04.23" "X1992.04.24"
 [848] "X1992.04.25" "X1992.04.26" "X1992.04.27" "X1992.04.28" "X1992.04.29" "X1992.04.30" "X1992.05.01"
 [855] "X1992.05.02" "X1992.05.03" "X1992.05.04" "X1992.05.05" "X1992.05.06" "X1992.05.07" "X1992.05.08"
 [862] "X1992.05.09" "X1992.05.10" "X1992.05.11" "X1992.05.12" "X1992.05.13" "X1992.05.14" "X1992.05.15"
 [869] "X1992.05.16" "X1992.05.17" "X1992.05.18" "X1992.05.19" "X1992.05.20" "X1992.05.21" "X1992.05.22"
 [876] "X1992.05.23" "X1992.05.24" "X1992.05.25" "X1992.05.26" "X1992.05.27" "X1992.05.28" "X1992.05.29"
 [883] "X1992.05.30" "X1992.05.31" "X1992.06.01" "X1992.06.02" "X1992.06.03" "X1992.06.04" "X1992.06.05"
 [890] "X1992.06.06" "X1992.06.07" "X1992.06.08" "X1992.06.09" "X1992.06.10" "X1992.06.11" "X1992.06.12"
 [897] "X1992.06.13" "X1992.06.14" "X1992.06.15" "X1992.06.16" "X1992.06.17" "X1992.06.18" "X1992.06.19"
 [904] "X1992.06.20" "X1992.06.21" "X1992.06.22" "X1992.06.23" "X1992.06.24" "X1992.06.25" "X1992.06.26"
 [911] "X1992.06.27" "X1992.06.28" "X1992.06.29" "X1992.06.30" "X1992.07.01" "X1992.07.02" "X1992.07.03"
 [918] "X1992.07.04" "X1992.07.05" "X1992.07.06" "X1992.07.07" "X1992.07.08" "X1992.07.09" "X1992.07.10"
 [925] "X1992.07.11" "X1992.07.12" "X1992.07.13" "X1992.07.14" "X1992.07.15" "X1992.07.16" "X1992.07.17"
 [932] "X1992.07.18" "X1992.07.19" "X1992.07.20" "X1992.07.21" "X1992.07.22" "X1992.07.23" "X1992.07.24"
 [939] "X1992.07.25" "X1992.07.26" "X1992.07.27" "X1992.07.28" "X1992.07.29" "X1992.07.30" "X1992.07.31"
 [946] "X1992.08.01" "X1992.08.02" "X1992.08.03" "X1992.08.04" "X1992.08.05" "X1992.08.06" "X1992.08.07"
 [953] "X1992.08.08" "X1992.08.09" "X1992.08.10" "X1992.08.11" "X1992.08.12" "X1992.08.13" "X1992.08.14"
 [960] "X1992.08.15" "X1992.08.16" "X1992.08.17" "X1992.08.18" "X1992.08.19" "X1992.08.20" "X1992.08.21"
 [967] "X1992.08.22" "X1992.08.23" "X1992.08.24" "X1992.08.25" "X1992.08.26" "X1992.08.27" "X1992.08.28"
 [974] "X1992.08.29" "X1992.08.30" "X1992.08.31" "X1992.09.01" "X1992.09.02" "X1992.09.03" "X1992.09.04"
 [981] "X1992.09.05" "X1992.09.06" "X1992.09.07" "X1992.09.08" "X1992.09.09" "X1992.09.10" "X1992.09.11"
 [988] "X1992.09.12" "X1992.09.13" "X1992.09.14" "X1992.09.15" "X1992.09.16" "X1992.09.17" "X1992.09.18"
 [995] "X1992.09.19" "X1992.09.20" "X1992.09.21" "X1992.09.22" "X1992.09.23" "X1992.09.24"
 [ reached getOption("max.print") -- omitted 8498 entries ]

Wide to long data

library(lubridate)

#wide to long
SEUS_temp <- melt(allvalues_wide_latlon, id.vars = c("lat", "lon"))
colnames(SEUS_temp) <- c("lat", "lon", "date", "sst")
SEUS_temp[,date.edit := as_date(substr(date, 2, 11))][,year := year(date.edit)][,sst_avg := mean(sst, na.rm = T), by = year]
ggplot(data = SEUS_temp, aes(x=year, y = sst_avg)) +
  geom_point() +
  geom_smooth(method = lm) +
  theme_classic()

Compare to SODA temp values

readRDS(here::here("Data","Spp_master","spp_master_ztemp.rds"))

readRDS(here::here("Data","Spp_master","spp_master_ztemp.rds"))
LS0tCnRpdGxlOiAiT0lTU0lUIE5PQUEgREFUQSIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKSGVyZSwgSSdtIHRyeWluZyB0byBnZXQgYSBiZXR0ZXIgc2Vuc2Ugb2YgaG93IHdlbGwgdGhlIFNPREEgdGVtcCBkYXRhIHJlcHJlc2VudHMgdGhlIFNFVVMuIFdlIGZvdW5kIGEgc2lnbmlmaWNhbnQgZGVjcmVhc2UgaW4gdGVtcGVyYXR1cmUgZHVyaW5nIHRoZSBzdHVkeSBwZXJpb2QsIHdoaWNoIGlzIG5vdCB3aGF0IHBlb3BsZSB3aG8gYXJlIGZhbWlsaWFyIHdpdGggY29uZGl0aW9ucyBpbiB0aGlzIGFyZWEgc2F5IGhhcyBoYXBwZW5lZC4gTW9zdCBsaWtlbHksIHRoZSBTT0RBIG1vZGVsIGRvZXNuJ3QgcGVyZm9ybSB3ZWxsIGluIHNoYWxsb3cgbmVhciBzaG9yZSBhcmVhcyB3aGVyZSBzdXJ2ZXkgaXMgcGVyZm9ybWVkIGluIHNvdXRoZWFzdCBVUy4gVGhlcmVmb3JlLCB3ZSBjb21wYXJlZCBwZXJmb3JtYW5jZSB0byBPSVNTVCB0ZW1wLiBPdmVyYWxsIHBhdHRlcm4gaXMgc2ltaWxhciwgYW5kIG5vdCByZXByZXNlbnRhdGl2ZSBvZiB3aGF0J3MgYWN0dWFsbHkgaGFwcGVuaW5nIGluIHRoaXMgc3VydmV5IGFyZWEuIFRoZXJlZm9yZSwgd2UgZGVjaWRlZCB0byBnbyB3aXRoIGluIHNpdHUgYnVveSB0ZW1wZXJhdHVyZSB2YWx1ZXMuIAoKYGBge3Igc2V0dXB9CmxpYnJhcnkobmNkZjQpCmxpYnJhcnkocmFzdGVyKQpsaWJyYXJ5KGRhdGEudGFibGUpCmBgYApTb3V0aGVhc3QgVVMgKDE5OTDigJMyMDE0LCBuID0gMjUpCgpZZWFycyBkb3dubG9hZGVkIG9uZSBieSBvbmUgb24gTWF5IDE1LCAyMDIwCmh0dHBzOi8vcHNsLm5vYWEuZ292L2NnaS1iaW4vRGF0YUFjY2Vzcy5wbD9EQl9kYXRhc2V0PU5PQUErSGlnaC1yZXNvbHV0aW9uK0JsZW5kZWQrQW5hbHlzaXMmREJfdmFyaWFibGU9U2VhK1N1cmZhY2UrVGVtcGVyYXR1cmUmREJfc3RhdGlzdGljPU1lYW4mREJfdGlkPTg2MTA3JkRCX2RpZD0xMzImREJfdmlkPTI0MjMKCmBgYHtyIGltcG9ydCBuZXRjZGZ9CiNhbGwgZmlsZSBuYW1lcwpmaWxlbmFtZXMgPC0gbGlzdC5maWxlcyhoZXJlOjpoZXJlKCJEYXRhIiwiVGVtcGVyYXR1cmUiLCJPSVNTVF9OT0FBX1RFTVBfTkVUQ0RGIiksIGZ1bGwubmFtZXMgPSBUKQoKT0lTU1Quc3N0IDwtIGxhcHBseShmaWxlbmFtZXMsIGJyaWNrKQoKT0lTU1Quc3N0CgoKI2N1cnJlbnRseSwgdGhpcyBpcyBpbiAwLTM2MCwgbmVlZHMgdG8gYmUgLTE4MC0xODAsIHNvIEkgaGF2ZSB0byBzdWJ0cmFjdCAzNjDLmgpPSVNTVC5zc3RfYWRqIDwtIGxhcHBseShPSVNTVC5zc3QsIHJvdGF0ZSkKCmBgYAoKRXh0cmFjdCBzcGVjaWZpYyBHUFMgdmFsdWVzIGZvciBTRSBVUyBzYW1wbGVkIGR1cmluZyBkdXJpbmcgc3R1ZHkgcGVyaW9kCmBgYHtyIFNFIFVTIHRlbXAgdmFsdWVzfQpsb2FkKGhlcmU6OmhlcmUoIkRhdGEiLCJTcHBfb2JzZXJ2YXRpb25zIiwidHJhd2xfY29vcmRpbmF0ZXNfMTk5Ni0yMDAwLlJEYXRhIikpCgpzZXVzX3BvaW50cyA8LSBncHNwb2ludHNvbmx5X3F1YXJ0ZXJkZWdyZWVyb3VuZFtyZWdpb24gPT0gIlNDRE5SX1NFVVMiXQoKY29vcmRpbmF0ZXMoc2V1c19wb2ludHMpID0gfiBsb24gKyBsYXQKCgp0ZW1wX3ZhbHVlc18xOTkwIDwtIGFzLmRhdGEudGFibGUoZXh0cmFjdChPSVNTVC5zc3RfYWRqW1sxXV0sIHNldXNfcG9pbnRzKSkKdGVtcF92YWx1ZXNfMTk5MSA8LSBhcy5kYXRhLnRhYmxlKGV4dHJhY3QoT0lTU1Quc3N0X2FkaltbMl1dLCBzZXVzX3BvaW50cykpCnRlbXBfdmFsdWVzXzE5OTIgPC0gYXMuZGF0YS50YWJsZShleHRyYWN0KE9JU1NULnNzdF9hZGpbWzNdXSwgc2V1c19wb2ludHMpKQp0ZW1wX3ZhbHVlc18xOTkzIDwtIGFzLmRhdGEudGFibGUoZXh0cmFjdChPSVNTVC5zc3RfYWRqW1s0XV0sIHNldXNfcG9pbnRzKSkKdGVtcF92YWx1ZXNfMTk5NCA8LSBhcy5kYXRhLnRhYmxlKGV4dHJhY3QoT0lTU1Quc3N0X2FkaltbNV1dLCBzZXVzX3BvaW50cykpCnRlbXBfdmFsdWVzXzE5OTUgPC0gYXMuZGF0YS50YWJsZShleHRyYWN0KE9JU1NULnNzdF9hZGpbWzZdXSwgc2V1c19wb2ludHMpKQp0ZW1wX3ZhbHVlc18xOTk2IDwtIGFzLmRhdGEudGFibGUoZXh0cmFjdChPSVNTVC5zc3RfYWRqW1s3XV0sIHNldXNfcG9pbnRzKSkKdGVtcF92YWx1ZXNfMTk5NyA8LSBhcy5kYXRhLnRhYmxlKGV4dHJhY3QoT0lTU1Quc3N0X2FkaltbOF1dLCBzZXVzX3BvaW50cykpCnRlbXBfdmFsdWVzXzE5OTggPC0gYXMuZGF0YS50YWJsZShleHRyYWN0KE9JU1NULnNzdF9hZGpbWzldXSwgc2V1c19wb2ludHMpKQp0ZW1wX3ZhbHVlc18xOTk5IDwtIGFzLmRhdGEudGFibGUoZXh0cmFjdChPSVNTVC5zc3RfYWRqW1sxMF1dLCBzZXVzX3BvaW50cykpCnRlbXBfdmFsdWVzXzIwMDAgPC0gYXMuZGF0YS50YWJsZShleHRyYWN0KE9JU1NULnNzdF9hZGpbWzExXV0sIHNldXNfcG9pbnRzKSkKdGVtcF92YWx1ZXNfMjAwMSA8LSBhcy5kYXRhLnRhYmxlKGV4dHJhY3QoT0lTU1Quc3N0X2FkaltbMTJdXSwgc2V1c19wb2ludHMpKQp0ZW1wX3ZhbHVlc18yMDAyIDwtIGFzLmRhdGEudGFibGUoZXh0cmFjdChPSVNTVC5zc3RfYWRqW1sxM11dLCBzZXVzX3BvaW50cykpCnRlbXBfdmFsdWVzXzIwMDMgPC0gYXMuZGF0YS50YWJsZShleHRyYWN0KE9JU1NULnNzdF9hZGpbWzE0XV0sIHNldXNfcG9pbnRzKSkKdGVtcF92YWx1ZXNfMjAwNCA8LSBhcy5kYXRhLnRhYmxlKGV4dHJhY3QoT0lTU1Quc3N0X2FkaltbMTVdXSwgc2V1c19wb2ludHMpKQp0ZW1wX3ZhbHVlc18yMDA1IDwtIGFzLmRhdGEudGFibGUoZXh0cmFjdChPSVNTVC5zc3RfYWRqW1sxNl1dLCBzZXVzX3BvaW50cykpCnRlbXBfdmFsdWVzXzIwMDYgPC0gYXMuZGF0YS50YWJsZShleHRyYWN0KE9JU1NULnNzdF9hZGpbWzE3XV0sIHNldXNfcG9pbnRzKSkKdGVtcF92YWx1ZXNfMjAwNyA8LSBhcy5kYXRhLnRhYmxlKGV4dHJhY3QoT0lTU1Quc3N0X2FkaltbMThdXSwgc2V1c19wb2ludHMpKQp0ZW1wX3ZhbHVlc18yMDA4IDwtIGFzLmRhdGEudGFibGUoZXh0cmFjdChPSVNTVC5zc3RfYWRqW1sxOV1dLCBzZXVzX3BvaW50cykpCnRlbXBfdmFsdWVzXzIwMDkgPC0gYXMuZGF0YS50YWJsZShleHRyYWN0KE9JU1NULnNzdF9hZGpbWzIwXV0sIHNldXNfcG9pbnRzKSkKdGVtcF92YWx1ZXNfMjAxMCA8LSBhcy5kYXRhLnRhYmxlKGV4dHJhY3QoT0lTU1Quc3N0X2FkaltbMjFdXSwgc2V1c19wb2ludHMpKQp0ZW1wX3ZhbHVlc18yMDExIDwtIGFzLmRhdGEudGFibGUoZXh0cmFjdChPSVNTVC5zc3RfYWRqW1syMl1dLCBzZXVzX3BvaW50cykpCnRlbXBfdmFsdWVzXzIwMTIgPC0gYXMuZGF0YS50YWJsZShleHRyYWN0KE9JU1NULnNzdF9hZGpbWzIzXV0sIHNldXNfcG9pbnRzKSkKdGVtcF92YWx1ZXNfMjAxMyA8LSBhcy5kYXRhLnRhYmxlKGV4dHJhY3QoT0lTU1Quc3N0X2FkaltbMjRdXSwgc2V1c19wb2ludHMpKQp0ZW1wX3ZhbHVlc18yMDE0IDwtIGFzLmRhdGEudGFibGUoZXh0cmFjdChPSVNTVC5zc3RfYWRqW1syNV1dLCBzZXVzX3BvaW50cykpCnRlbXBfdmFsdWVzXzIwMTUgPC0gYXMuZGF0YS50YWJsZShleHRyYWN0KE9JU1NULnNzdF9hZGpbWzI2XV0sIHNldXNfcG9pbnRzKSkKCnBsb3QoT0lTU1Quc3N0X2FkaltbMV1dW1sxXV0pCnBvaW50cyhzZXVzX3BvaW50cykKCm9iamVjdF9uYW1lcyA8LSBsaXN0KHRlbXBfdmFsdWVzXzE5OTAsIHRlbXBfdmFsdWVzXzE5OTEsIHRlbXBfdmFsdWVzXzE5OTIsIHRlbXBfdmFsdWVzXzE5OTMsIHRlbXBfdmFsdWVzXzE5OTQsIHRlbXBfdmFsdWVzXzE5OTUsIHRlbXBfdmFsdWVzXzE5OTYsIHRlbXBfdmFsdWVzXzE5OTcsIHRlbXBfdmFsdWVzXzE5OTgsIHRlbXBfdmFsdWVzXzE5OTksIHRlbXBfdmFsdWVzXzIwMDAsIHRlbXBfdmFsdWVzXzIwMDEsIHRlbXBfdmFsdWVzXzIwMDIsIHRlbXBfdmFsdWVzXzIwMDMsIHRlbXBfdmFsdWVzXzIwMDQsIHRlbXBfdmFsdWVzXzIwMDUsIHRlbXBfdmFsdWVzXzIwMDYsIHRlbXBfdmFsdWVzXzIwMDcsIHRlbXBfdmFsdWVzXzIwMDgsIHRlbXBfdmFsdWVzXzIwMDksIHRlbXBfdmFsdWVzXzIwMTAsIHRlbXBfdmFsdWVzXzIwMTEsIHRlbXBfdmFsdWVzXzIwMTIsIHRlbXBfdmFsdWVzXzIwMTMsIHRlbXBfdmFsdWVzXzIwMTQsIHRlbXBfdmFsdWVzXzIwMTUpCgphbGx2YWx1ZXNfd2lkZSA8LSBkby5jYWxsKCJjYmluZCIsIG9iamVjdF9uYW1lcykKCmFsbHZhbHVlc193aWRlX2xhdGxvbiA8LSBjYmluZChzZXVzX3BvaW50c0Bjb29yZHMsIGFsbHZhbHVlc193aWRlKQpjb2xuYW1lcyhhbGx2YWx1ZXNfd2lkZV9sYXRsb24pCgpgYGAKCldpZGUgdG8gbG9uZyBkYXRhCmBgYHtyIHdpZGUgdG8gbG9uZ30KbGlicmFyeShsdWJyaWRhdGUpCgojd2lkZSB0byBsb25nClNFVVNfdGVtcCA8LSBtZWx0KGFsbHZhbHVlc193aWRlX2xhdGxvbiwgaWQudmFycyA9IGMoImxhdCIsICJsb24iKSkKY29sbmFtZXMoU0VVU190ZW1wKSA8LSBjKCJsYXQiLCAibG9uIiwgImRhdGUiLCAic3N0IikKU0VVU190ZW1wWyxkYXRlLmVkaXQgOj0gYXNfZGF0ZShzdWJzdHIoZGF0ZSwgMiwgMTEpKV1bLHllYXIgOj0geWVhcihkYXRlLmVkaXQpXVssc3N0X2F2ZyA6PSBtZWFuKHNzdCwgbmEucm0gPSBUKSwgYnkgPSB5ZWFyXQoKCgpgYGAKCgpgYGB7ciB0ZW1wIG92ZXIgdGltZX0KZ2dwbG90KGRhdGEgPSBTRVVTX3RlbXAsIGFlcyh4PXllYXIsIHkgPSBzc3RfYXZnKSkgKwogIGdlb21fcG9pbnQoKSArCiAgZ2VvbV9zbW9vdGgobWV0aG9kID0gbG0pICsKICB0aGVtZV9jbGFzc2ljKCkKYGBgCgpDb21wYXJlIHRvIFNPREEgdGVtcCB2YWx1ZXMKCmBgYHtyIHNvZGEgdmVyc3VzIE9JU1NUfQpzcHBfbWFzdGVyX3p0ZW1wIDwtIHJlYWRSRFMoaGVyZTo6aGVyZSgiRGF0YSIsIlNwcF9tYXN0ZXIiLCJzcHBfbWFzdGVyX3p0ZW1wLnJkcyIpKQoKZ2dwbG90KGRhdGEgPSBTRVVTX3RlbXAsIGFlcyh4PXllYXIsIHkgPSBzc3RfYXZnKSkgKwogIGdlb21fcG9pbnQoY29sb3IgPSAicmVkIikgKwogIGdlb21fc21vb3RoKG1ldGhvZCA9IGxtLCBjb2xvciA9ICJyZWQiKSArCiAgZ2VvbV9wb2ludChkYXRhID0gc3BwX21hc3Rlcl96dGVtcFtyZWcgPT0gInNhIl0sIGFlcyh4PXllYXIsIHkgPSBtZWFuX3NzdF90ZW1wKSwgY29sb3IgPSAiYmx1ZSIpICsKICBnZW9tX3Ntb290aChkYXRhID0gc3BwX21hc3Rlcl96dGVtcFtyZWcgPT0gInNhIl0sIGFlcyh4PXllYXIsIHkgPSBtZWFuX3NzdF90ZW1wKSwgbWV0aG9kID0gbG0sIGNvbG9yID0gImJsdWUiKSArCiAgdGhlbWVfY2xhc3NpYygpCmBgYAoK